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Abstract 

We present measurements of the microwave surface resistance Rs and the penetra- 
tion depth A of YiBa2Cu307„5 crystals. At low T, A(T) obeys a polynomial be- 
havior, while -Rs displays a characteristic non-monotonic T— dependence. A detailed 
comparison of the experimental data is made to a model of d-wave superconductivity 
which includes both elastic and inelastic scattering. While the model reproduces the 
general features of the experimental data, three aspects of the parameters needed 
are worth noting. The elastic scattering rate required to fit the data is much smaller 
than measured from the normal state, the scattering phase shifts have to be close to 
7r/2 and a strong coupling value of the gap parameter 2A(0)/kTc ~ 6 is needed. On 
the experimental side the uncertainties regarding the material parameters A(0) and 
-Rs,res(0) further complicate a quantitative comparison. For one sample, -Rs,res(0) 
agrees with the intrinsic value which results from the d-wave model. 



Microwave measurements of the surface impedance = Rs + iXs of super- 
conductors are in principle capable of yielding a wealth of precise information 
regarding the superconducting state, such as the gap parameter, quasiparticle 
density and nature of scattering. In low Tc superconductors the BCS theory 
provides a remarkably accurate description of experimental data for Rg and 
Xf, over several orders of magnitude variation, including detailed effects of 
impurity scattering [1]. 
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Table 1 

Some sample properties: Tc is observed by us in small microwave fields. AT^ is 
defined as the temperature interval between 10% and 90% of Rs{Tc). S is the oxygen 
vacancy. 



Sample 




Ar, 


-Rs(-iK) 


5 


A 


91.0K 


l.OK 




0.06 


B 


91.5K 


1.5K 


150 


0.11 



Recently, experiments which directly explore the order parameter symmetry 
suggest a dr^2_y2 order parameter [2-4] for high Tc superconductors, although 
some experiments which suggest a s-wave order parameter also exist [5,6]. It 
is therefore useful to ask to what extent a d-wave model of superconductivity 
can describe the measured surface impedance of the cuprate superconductors, 
particularly YiBa2Cu307_5. Some work has been already initiated in this re- 
gard [7]. Here in this paper, we present detailed results on the microwave 
(10 GHz) surface impedance of YBa2Cu307_5 crystals. We also compare the 
complete temperature dependence to numerical calculations based upon semi- 
microscopic models of (i-wavc and s-wave superconductivity, including elastic 
as well as inelastic scattering effects. 

The measurements were carried out in a specially designed, high sensitivity Nb 
cavity. The method of measuring the surface impedance of superconductors at 
elevated temperatures using a "hot finger" cavity method was first introduced 
by one of the authors in reference [8], and the principle has been used in a 
variety of systems reported in the literature. In the present setup, the Nb 
cavity is maintained either at 4.2 K or below 2K. The typical background 
Qi, of the cavity can be as high as 10^. The surface resistance is measured 
from the temperature dependent Q using Rs{T) = T [Q~^{T) — Qj^^(T)] and 
the penetration depth using AA(T) = ( [f(T) — fb(T)]. The geometric factors 
are determined by the cavity mode, sample location and the sample size. All 
measurements were done in the TEqh mode with the sample at the midpoint 
of the cavity axis, where the microwave magnetic fields have a maximum and 
the microwave electric fields are zero. The method enables measurement of 
small crystals and thin film samples. 

The crystals were grown in a Z1CO2/Y crucible using highly pure Y2O3, BaCOa 
and CtiO powders. Crystal growth took place while slowly cooling the melt 
at a rate of 0.40°C/h in the temperature range 970°C to 904°C and in an 
atmosphere of lOOmbar O2. After the growth the crystals were annealed in 
flowing oxygen in the temperature range 600°C to 400°C during 600 h. 

We start with an analysis of the normal state properties in the hope to fix 
some of the normal state parameters which are required for a calculation of the 
conductivity in the superconducting state. Assuming local electrodynamics, 
i.e. skin depth limited, the normal state surface resistance is given by 



2 




200 I ' < ' < ' < ' < ' < ' < ' < ' < ' < ' 1 

100 120 140 160 180 200 220 240 260 280 300 

T[K] 

Fig. 1. Normal state surface resistance RsiT) for the YBa2Cu307_5 crystals with 
fits to eq. (2). The fit parameters are given in table 2. 

Table 2 

Experimental normal state scattering parameters 



Sample 


Po 


7 


Tel 


a 






[nm/K] 


[sec--^] 


[sec-i K-i] 


A 


5.310"^ 


9.110"^ 


1.110^3 


1.910" 


B 


1.410-6 


2.410-8 


2.910^3 


5.010" 



i?n = where pn = o'cT^ = /^0'^(0)^2r = (ne^r/m)-^. If the classical 

skin depth limit applies, then the microwave resistivity pn should be the same 
as the dc resistivity, which is known to be linear. If Pn = Po + iT, then 

i?n = v/^po(po + 72^)/2 (1) 

Figure 1 shows a comparison of the experimental data between 100 K and 
200 K to the above equation. The data clearly has a sub-linear T dependence, 
and the fit to eq. (1) is extremely good. From the fit the scattering rate F can 
be obtained as 

F = Fei + Finei(T) = Fei + aT (2) 

where Fei = po/(2po'^(0)^), and a = 7/(2/xoAq). Table 2 gives these parameters 
for the studied samples under the assumption A(0) = 1400 A. 

Figure 2 displays the low temperature behavior of Rg for the two samples of 
YiBa2Cu307_5. Sample B shows a characteristic peak in first reported by 
[9]. This is not present in sample A, whose behavior is instead closer to that 
of thin films. 

It is evident that the A vs. T data displayed in figure 4 do not show the 
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Fig. 2. Low temperature behavior of Rg for the YBa2Cu30i„5 samples A (•) and 
B (A). Theoretical results are plotted for the parameters given in table 3. 

Table 3 

Theoretical parameter with fi{t) = e'^^*~^^ — e~'^ 



Sample 


symbol 


Pel [meV] 


Find [meV] 


2Ao(0)/kre 


A 


dashes 


0.25 


10.86 


6.4 




dots 


0.20 


10.91*3 


6.0 


B 


dash dot 


0.05 


24.15 


6.4 




line 


0.05 


24.15*3 


7.0 



exponential dependence AA(T) oc exp(— T/Tc) expected for an isotropic s- 
wave superconductor. Instead the data clearly have a polynomial temperature 
dependence with a leading linear term [10,11]. This, together with the nonex- 
ponential decrease of Rs at very low temperatures could be taken as indication 
for d-wave pairing. 

An important prediction for any superconducting state with nodes in the gap 
is the presence of a finite residual conductivity uoo due to elastic scattering 
[12,13]. Its value depends on the particular pair state and, if determined ex- 
perimentally, could help to identify the type of pairing present. For the d-wave 
state 

A(r, kp) = Ao(T) cos(2</.) (3) 

commonly studied for systems with cylindrical Fermi surfaces of circular cross 
section one has [12,13] aoo = ne^/m7rAo at T = OK. 

This can be related to o'n(Tc) by o-oo/'^ii(^c) = 2r(Tc)/7rAo. Using the values 
2r(Tc) = 33 meV, and 2Ao/kTc = 4.3 we get o'oo/o'n(^c) ~ 1- Thus the residual 
conductivity is comparable to the normal state conductivity at T^, which is a 
surprisingly large value. Note that in BCS the conductivity (7i — > as T — > 0. 
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The relation to Rg is obtained from the hmiting result when o"2 ^ o"i, whence 

R^(T ^ 0) = cj2/i2^(0)Voo/2. This can also be written as Rs/Rn = 0.5((Ti/an)/(cT2/(Jn)^/^. 

Since the above relationship shows that ai/a^ ~ 1, the reduction in Rg from 

its value at Tc is entirely due to the reduction in A from 5n i- e- to the superfiuid 

response. The intrinsic residual surface resistance i?s,res(0) which results from 

(Too is much smaller than the value measured for sample B but is compatible 

with the low temperature data for sample A. 

Detailed numerical calculations based on the microscopic models for s-wave 
and d-wave superconductivity [14-16] were carried out. 

With suitable choices for the Eliashberg functions it is possible to fit the mea- 
sured surface resistances within the framework of isotropic strong coupling 
theory over a wide range of temperature from Tc down to 0.4 Tc [16], high- 
lighting the importance of inelastic scattering processes. This is consistent 
with earlier measurements over the same temperature range [17]. However at 
lower temperatures the isotropic gap should make its presence known in the 
form of an exp(— A/kT) dependence, which is clearly not observed in either 
i?s or AA. 

For this reason we focus in this paper on models of d-wave superconductivity, 
including inelastic scattering. Here we list the main ingredients of the model 
- details of the model are described elsewhere [16]. A weak-couphng pairing 

— * 

interaction was suitably chosen to give the dx2-y2 order parameter A(T, /cp), 
see eq. (3). The self-consistency equation for the order parameter was solved 
to give the temperature dependence of the amplitude Ao(T), and which is 
found to be very similar to that for an isotropic order parameter except that 
2Ao(0)/kTc = 4.29 rather than 3.52. 

Elastic scattering is parametrized by a normal state scattering rate Fd and 
a phase shift which can take any value between (Born approximation) 
and 7t/2 (Unitary limit). Inelastic scattering in the superconducting state is 
parametrized in the form 

Finei(r) = aTJit) (4) 

with some function / of the reduced temperature t — T/Tf.. 

The inputs to the calculation of the surface impedance Zs = (2i/cTs)^/^ are 
then Fgi, ct, and f{t). In order to fit the steep drop of i?s below Tc we have 
found it necessary to treat 2Ao(0)/kTc as a variable parameter. 

The surface impedance is not as sensitive to the choice of model parameters 
as is the real part of the conductivity ai, which is also of some intrinsic inter- 
est. Unfortunately, the peak height of the experimentally determined (Ti(T, a;) 
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Fig. 3. Real part ai of the measured conductivity and theoretical behavior (lines). 
The parameters are given in table 3. 

depends strongly on the choice of the zero temperature penetration depth 
A(0) while the low temperature behavior of ai can be changed substantially 
by subtracting from the measured surface resistance a residual loss R° [9]. 
The choice of R° is limited by the consideration that o"i(T = 0,uj) should 
neither be negative nor should it exceed (Ti(Tc,a;) by a wide margin. With 
these restrictions we find it!° = 0.15 mQ for sample B while for sample A R° 
is neghgible. 

The peak heights in o"i are reduced when larger values for A(0) are selected. We 
have chosen A(0) such that the experimental Xg is linear over as large a temper- 
ature range below Tc as possible. This gives A(0) = 1800 A and a^^{T,uj) — 
30 [mQcm]-^ for sample A and A(0) = 2000 A and a'^^{T,uj) = 42 [mQcm]-^ 
for sample B. 



Experimental results for cti are shown in figure 3. In the presence of purely 
inelastic scattering a peak in ai should occur near a;/2r(T) = 1. This peak 
is thus expected to shift to higher temperatures as the frequency is increased 
and to lower temperatures when the overall magnitude a (eq. (2) and (4)) of 
the scattering rate is increased. For f(t) = fi{t) = e^^*"^-* — e~'^ which has 
been suggested by Bonn et al. [9] and which is close to the result from the 
Nested Fermi Liquid model [18,16], the peak in cti occurs near 10 K with the 
peak height greatly exceeding the experimental value. This observation, taken 
together with the normal state data, shows that elastic scattering needs to 
be taken into account. When Fd given in table 2 is used as input parameter, 
the peak in cxi is greatly reduced. In the Born approximation cxi rises very 
steeply from its limiting value (Too so that the peak is still located at too 
low a temperature. In the unitary limit a peak is barely observable. In order 
to reproduce the experimental results, Fd has to be chosen much smaller 
than the analysis of normal state data would suggest. Figure 3 contains a 
reasonably close fit to the Ui data. The fit cannot be improved by varying the 
phase shift 5n- Reducing 5n from O.Stt shifts some of the weight of the peak 
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Fig. 4. Low temperature data and theory of A(r). Theoretical curves are obtained 
with the parameters of table 3. 

shown in figure 3 to the temperature at which the peak occurs in the Born 
approximation. At around 6-^ = 0.357r, cri(T) acquires a distinct double peak 
structure not compatible with the data. We conclude that the phase shift must 
be close to the unitary hmit 0.457r < 5n < O.Stt. 

The only way to improve the fit would be to choose different temperature 
dependencies for Finei with f{t) decreasing faster then = e^*-*^^'' — for 
sample B and more slowly for sample A. Even though the two samples differ 
in oxygen contents (table 1) and by a factor of four in thickness, it does not 
seem plausible that intrinsic scattering events in the two samples should differ 
significantly in their temperature dependencies. A more likely source for this 
discrepancy is a temperature dependence of the "residual" surface resistance 
Rs,ies{t)- Models have been put forward to explain i?s,res(^) by weak links and 
relate them to po [19]. 

The theoretical surface resistance for the same model parameters as those used 
to calculate ai is shown in figure 2. In the case of sample B the phenomenolog- 
ical residual resistance i?s,res(0) = 0.15 mil has been added. Note the intrinsic 
residual surface resistance in the case of sample A. To fit the data near Tc 
we had to increase 2Ao(0)/kTc to 6.4. This is substantially larger than the 
weak coupling value 4.29, which has important implications for conclusions 
regarding fluctuations. The overall fit to the data is very good, showing the 
same small discrepancies already apparent in figure 3. 

The case for (i-wave pairing would be strengthened considerably if we could fit 
the shift in penetration depth equally well using exactly the same model. Ex- 
perimental and theoretical results are compared in figure 4. Clearly, the agree- 
ment is less than satisfactory. Note that the calculated AA is by no means lin- 
ear, although a polynomial fit in a limited temperature range can certainly be 
found. The agreement could be much improved by choosing different tempera- 
ture dependences for Finei- A good fit for sample A is obtained with /2(t) = t^, 
see figure 4. Decreasing 2Ao(0)/kTc substantially would also improve agree- 
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Fig. 5. Superfluid density of sample A. Dashed and doted lines are obtained with 
the parameters of table 3. 

ment in the case of AA but would lead to serious discrepancies in the case of 

In figure 5 the penetration depth for sample A is plotted over the whole 
temperature range in the form of the superfluid density (A(0)/A(T))^. For 
such a plot, it is necessary to assume a value of A(0) and this figure shows the 
variation of (A(0)/A(T))^ resulting from different choices of A(0). By varying 
A(0) one can actually change the sign of the curvature of A^(T)) near Tc. 
It is obvious that near Tc, since a straight line fits the data pretty well, the 
behavior is described quite well by a (1 — T/Tc)^/^ dependence, which suggests 
a mean-field behavior of the order parameter. 

Our (i-wave calculations give a positive curvature for 1/A^(T) near T^, indi- 
cating a behavior which is slower than mean-field, i.e. A(T) (1 —t)", where 
u > 1/2. This is possibly an artifact of the calculation because the total scat- 
tering rate (eq. (2)) at Tc is such that Tc should be substantially suppressed 
according to weak coupling theory. This would lead to noticeably different Tc's 
for the two samples. Since the Tc's are practically the same, we did not include 
scattering in the selfconsistency equation. Strong inelastic scattering proba- 
bly has to be treated within the framework of an anisotropic strong coupling 
theory, which could also solve the problem of the large value for 2Ao(0)/kTc 
we had to assume. 

In spite of the remaining discrepancies, some of which may be due to con- 
tributions from the c-axis conductivity which has not been included in the 
calculations, d-wave pairing seems to provide an adequate model for under- 
standing features seen in YBa2Cu307_5 crystals. The main features of the data 
appear to be reproduced, although a detailed microscopic justification of the 
needed parameters is not yet available. We should remark that although we 
have considered an explicit d-wave model, the essential feature is that of nodes 
in the gap leading to low lying quasiparticle excitations at all temperatures. 
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